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The white dwarf luminosity function, which provides information about their coohng, has 

>- been measured with high precision in the past few years. Simulations that include well 

^ known Standard Model physics give a good fit to the data. This leaves little room for new 

^ physics and makes these astrophysical objects a good laboratory for testing models beyond 

(T^ the Standard Model. It has already been suggested that white dwarfs might provide some 

^ evidence for the existence of axions. In this work we study the constraints that the white 
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dwarf luminosity function puts on physics beyond the Standard Model involving new light 
particles (fermions or bosons) that can be pair-produced in a white dwarf and then escape 
to contribute to its cooling. We show, in particular, that we can severely constrain the 
parameter space of models with dark forces and light hidden sectors (lighter than a few 
keV). The bounds we find are more competitive than those from current lab searches and 
those expected from most future searches. 
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1. Introduction 



White dwarfs (WDs) are simple astrophysical objects wliose cooling law is well understood. 
This fact makes them a good laboratory for testing new models of particle physics. Many 
such models predict the existence of light bosons or light fermions that interact very weakly 
with regular matter. If these new particles are produced in a WD, they will typically escape 
and accelerate the cooling of the star. Thus, determining the cooling law from astrophysical 
observations can be translated into constraints on particle physics beyond the Standard 
Model (BSM) 

We first give a brief review of WD cooling. Formally, the cooling evolution of white 
dwarfs can be written as [2]: 

L, + L. + L., = -y^ C^dm-j^ T[g^)^^^^^dm+il. + e,)— + ex, (1.1) 

where and Li, represent the photon and neutrino luminosities (energy per unit time). 
The first term on the r.h.s. is the well known contribution of the heat capacity of the star 
to the total luminosity, the second one represents the contribution of the change of volume. 
It is in general small since only the thermal part of the electronic pressure, the ideal part 
of the ions and the Coulomb terms other than the Madelung term contribute. The third 
term represents the contribution of the latent heat and gravitational readjustement of the 
white dwarf to the total luminosity at freezing. Finally, and represent any extra 
energy sink or source of energy respectively. For many applications, this equation can be 
easily evaluated assuming an isothermal, almost completely degenerate core containing the 
bulk of mass, surrounded by a thin, nondegenerate envelope]^ 

The evolution of white dwarfs can be tested through the luminosity function (LF), 
n{l), which is defined as the number of white dwarfs of a given luminosity or bolometric 



^The isothermal approximation is not valid when neutrinos are dominant, however the results are still 
reasonably good and provide a reasonable guess of the luminosity. The ions do not follow the ideal gas law but 
the equation of state of a Coulomb plasma — for instance, in the region of interest the specific heat approaches 
the Dulong-Petit law — and crystallizes at low temperatures, around bolometric magnitude 12 — 13, depending 
on the mass of the star. 
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magnitud^ per unit of magnitude interval and unit volume: 

/•Ms 

n{l) = / $(M) ^(r)reooi(/, M) dM (1.3) 

J Mi 

where 

r = TG-teool(/,M)-tpg(M). (1.4) 

Tg is the age of the Galaxy, / = — log(L/L0), M is the mass of the parent star (for 
convenience all white dwarfs are labeled with the mass of the main sequence progenitor), 
tcooi is the cooling time down to luminosity /, TcooI = dt/dM\,o\ is the characteristic cooling 
time, Ms is the maximum mass of a main sequence star able to produce a white dwarf, 
and Mi is the minimum mass of the main sequence stars able to produce a white dwarf 
of luminosity /, and tpg is the lifetime of the progenitor of the white dwarf. $(M) is the 
initial mass function, i.e. the number of main sequence stars of mass M that are born per 
unit mass, and \l/(t) is the star formation rate, i.e. the mass per unit time and volume 
converted into stars. So the product $(M)\E'(r) is the number of main sequence stars 
that were born at the right moment to produce a white dwarf of luminosiy / now. Since 
the total density of white dwarfs is not well known, the computed luminosity function is 
usually normalized to the bin with the smallest error bar, traditionally the one with / = 3, 
in order to compare theory with observations. 



The star formation rate is not known but fortunately the bright part of Eq. (1.3) 
satisfies [3]: 

n{l) oc (tcooi) I $(M)^(r) dM 



If \E' is a well behaved function and Tq is large enough, the lower limit of the integral is 
not sensitive to the luminosity, and its value is absorbed by the normalization procedure 
in such a way that the shape of the luminosity function only depends on the averaged 
characteristic cooling time of white dwarfs. 

It is important to realize that white dwarfs are divided into two broad categories, DA 
and non-DA. The DA white dwarfs exhibit hydrogen lines in their spectra caused by the 
presence of an external layer made of almost pure H. This hydrogen layer is absent in the 
case of non-DAs and, consequently, their spectra is free of the H spectral features. The 
main result is that the DAs cool down more slowly than the non-DAs |7j. 



^Thc bolometric magnitude and the luminosity are related through 

Mboi = -2.5 logio(i/i0) + 4.74. (1.2) 
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Fig. 1: Luminosity function of white dwarfs. Red |4| and blue [5] points represent tlie 
luminosity function of all white dwarfs (DA and non-DA families). Magenta points [6] 
represent the luminosity function of the DA white dwarfs alone. Both distributions have 
been normalized around Mboi = 13, see text. The dotted line represents the luminosity 
function obtained assuming Mestel's approximation. The continuous lines correspond to 
full simulations assuming a constant star formation rate and an age of the Galaxy of 13 
Gyr for the DA family (black line) and all, DA and non-DA, white dwarfs (blue line). 

The observed LF is shown in Fig. [T] for three different datasets. Note that moving 
from left to right along the horizontal axis we go from high luminosity (hot, young WDs) 
to low luminosity (cold, old WDs). The Harris et al. |4] and the Krzesinski et al. [5j are 
representative of all, DAs and non-DAS, white dwarfs. The Harris et al. LF has been 
constructed using the reduced proper motion method which is accurate for cold WDs with 
^boi ^ 6 but not appropriate for hot WDs with Mboi ^ 6, which vave been removed from 
the sample. The Krzesinski et al. LF on the other hand has been built employing the 
UV-excess technique which is accurate for hot WDs with Mboi ^ 7 but inappropriate for 
the colder ones. Since the datasets overlap, assuming continuity, it is possible to construct 
a LF that extends from Mboi ~ 1-5 to Mboi ~ 16, although the cool end is affected by 
severe selection effects. The DeGennaro et al. [6] sample was also obtained with the proper 
motion technique, for which reason the hot end is not reliable and has been removed. 
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Since DAs and non-DAs confuse at low temperatures, the corresponding points have been 
removed from Figure [l] See Isern et al. [s] for a detailed discussion. Since ultimately only 
the slope of the LF is of interest and the total density of WDs is quite uncertain, it is 
usually more convenient to normalize the LF with respect to one of its values, which is 
commonly chosen around log(L/L0) = —3. 

If the cooling were due only to photons and one assumes that Mestel's approximation [o] 
holds {i.e. ions behave like an ideal gas and the opacity of the radiative envelope follows 
Kramer's law), then the LF would be a straight line on this logarithmic plot, which already 
provides a reasonable fit to the data. Note, however, that the data show a dip for values 
of Mboi around 6 — 7. That is where the neutrinos enter the game: for the hotter WDs (to 
the left in Fig. [T]), neutrino emission becomes more important than photon cooling. When 
neutrinos are included and the cooling is followed with a full stellar evolution code the 
agreement becomes impressive (see the continuous lines of Fig. [T]). This agreement can be 
used to bound the inclusion of new sources or sinks of energy Isl. 



2. Cooling mechanisms 

In this section we review the various cooling mechanisms for WDs. The aim is to provide 
the reader with a simple understanding of what mechanism dominates in what regime. 

2.1. Photons 

In WDs the thermal energy is mostly stored in the nuclei which form, to a good approxima- 
tion, a classical Boltzmann gas. Taking into account the thermal conductance of the surface 
layers, one can relate the rate of energy loss at the surface to the internal temperature. 
Using Mestel's approximation |9] one finds 

= 3.29 X 10"^ T//^ erg g"^ s"\ (2.1) 

where e-y is the energy-loss rate per unit mass and T-j = j^-^. This constitutes the main 
cooling for cold (Mboi ^ 7) WDs. Realistic models indicate that oc T^, where /3 ~ 7/2, 
but varies slightly with temperature, chemical composition and mass of the white dwarf. 
We show in Fig. |2] what this energy loss as a function of the core temperature looks like 
for a realistic model as opposed to Mestel's model. 
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Fig. 2: Luminosity (bolometric magnitude) versus core temperature for a realistic model 
(continuous line) and for Mestel's model (dashed line). The photon luminosity ~ 
e-yMwD, with Mwd the WD mass, is related to Mboi via Eq. (1.2). 



2.2. Light bosons vs light fermions in white dwarfs 

Additional light bosons and light fermions that interact very weakly can also contribute to 
the cooling of WDs, but their dominant production mechanisms are usually different. First, 



consider the DFSZ axion 10 as an example of a light boson. It would be mainly produced 
by the bremsstrahlung process e + (Z, A) ^ e + (Z, A) + a as shown in Fig. Is] (a) 



Raffelt 

gave an intuitive argument |11| to understand how the related energy emission rate depends 
on the temperature. It goes as follows: The relevant interaction term in the Lagrangian 
is igae'j^e, where g = me//pQ, with me the electron mass, /pq > 10^ GeV the Peccei- 
Quinn scale, a the axion field and e the electron field. The axion emission by an electron 
is analogous to the emission of a photon but, due to the presence of the 75, there is an 
extra electron spin- flip in the amplitude. Whereas the usual photon bremsstrahlung cross 
section is proportional to E~^, the axionic analogue is proportional to Ea due to the extra 
power from the spin-flip nature of the process. For the energy emission rate, we have to 
multiply the cross section by another factor of Ea, which makes it proportional to E^. We 
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Fig. 3: Processes for the production of new light scalars and fermions in WDs. The upper 
two diagrams represent bremsstraiilung processes, wliile tlie lower two diagrams show the 
plasmon decay. The latter is the main production mechanism in WDs when the new light 
particles are produced in pairs. The red dot in these diagrams represents an effective 
interaction given by a dimension 6 operator. The grey blob in the lower two diagrams 
represents the effects of the stellar medium, as discussed in the text. 

still have to do the phase space integrals for the initial and final state electrons. Because 
electrons are degenerate in WDs, these integrals contribute a factor of T/Ep each, with 
Ep the electron Fermi energy. At the end, the rate is proportional to E'^{T/EpY oc T^, 
given that axion energies will be of the order of the temperature T. 

Next, consider the bremsstrahlung of two fermions ■?/', as depicted in Fig. [3] (b). The 
intuitive reasoning is analogous to what we just described for axions, but the difference 
is that the fermions from the electron line are produced in pairs as opposed to the single 
axion. That adds an extra factor of the energy (~ T) in the cross section and an extra 
phase space integral. As a result, we have two more powers of T in the final rate, which 
is therefore proportional to T^. 

When the calculations are done carefully one gets the following results for the energy- 
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loss rates per unit mass in the two cases [T] 



1i 

,brem _ j Qg ^ y4 \^ F„ erg g"^ S"^ 



ej- = 1.34 X 10-^ (^)'Tr^ E^^f 

Here, = 10^^ with (7 the couphng of the axion to electrons defined abovej^ is the 
dimensionful coupling for the four-fermion interaction denoted by a red dot in Fig. [3| to 
be compared to the familiar Fermi constant, Gp = 1-166 x 10^^ GeV^^; is the effective 
coupling constant analogous to the effective neutral-current vector coupling constant Cy = 
0.964; Xj is the mass fraction of the element j, with nuclear charge Zj and atomic mass 
number Aj, and the sum runs over the species of nuclei present in the WD. Fa and are 
factors that take into account the effect of screening for Coulomb scattering in a plasma. 
In WDs Fa and F^ are of order one to a good approximation. 



It is clear from expressions (2.1) and (2.2) why the bremsstrahlung process for the 
neutrinos, where C^G^ = CyGp, is completely irrelevant in WDs, with internal temperature 
of the order of 10^ K. First, the numerical coefficient in e^""^™ is suppressed by four orders of 
magnitude compared to photons (and to axions if we take of order one). Second, it has 
a steeper dependence on the temperature, which makes it less and less relevant as we go to 
lower temperatures (see Fig. |4|. Unless we have a model in which is significantly bigger 
than Gp, this contribution is negligible. In fact, the dominant production mechanism of a 
pair of light fermions in WDs is not bremsstrahlung but is given by the so-called plasmon 
process [13,|, which is depicted in Fig. |3] (c). That is what we describe next. 

In vacuum the photon is massless and can not decay into a pair of massive particles, 
no matter how light they are. But in a medium, as in the interior of a star, the photon 
dispersion relations are modified and this allows such a decay. What happens is that the 
photon also acquires a longitudinal polarization and is promoted to the so-called plasmon. 
One would be tempted to say that the photon becomes massive, but such a statement is 
strictly speaking incorrect. A better way to think about the plasmon decay, without ever 
referring to the mass of the photon, is the following: the propagation of an electromagnetic 
excitation (the plasmon) in the plasma is accompanied by an organized oscillation of the 
electrons, that in turn serve as a source for emitting a pair of light particles. Figs. IS] (c,d) 



^For a26 of order one, axion cooling becomes comparable to photon cooling and one gets a better fit to 
the LF [1^. This fact can be taken as tentative evidence for the existence of axions, and it explains the 
choice of the power of 26 in the definition of a26- 
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Photon 

Plasmon (iptp) 

— Bremsstrahlung (a) 

— Bremsstrahlung {ipip) 



Fig. 4: Comparison of energy losses in WDs. In this plot, T-j = jqt-k, "we have set = 0.3, 
C^G^ = CyGp, and we have assumed a WD composed of an equal mixture of ^^C and ^^O. 
For two fermions in the final state, the plasmon contribution dominates at high internal 
temperature, Ty > 5, while the bremsstrahlung contribution is completely negligible. 



are then understood as follows: the grey blob represents the medium response to the elec- 
tromagnetic excitation; we can think of the black line outlining the blob as a loop of 
electrons, with the red dot denoting an effective interaction with the pair of light parti- 
cles, that can be either fermions or bosons. This is a schematic description. The reader 
interested in more details is referred to the pedagogic treatment in chapter 6 of Ref. [l]. 



The calculation of the plasmon decay 13 ,14 is quite involved, due to the effects of the 



medium, and cannot be performed analytically. However, a good approximation, in the 



case of neutrinos as the products of the decay, was given in Ref. 15 . The result applies to 



a wide range of stellar temperatures and densities. Restricting ourselves to WDs, we can 
write it as 

2 



plasmon 



.40 X 10^^ X'l'e-\fT + h) erg g-^ s'^ 



(2.3) 



where numerically, to a good approximation 



A = 1.69 X 10"^ T- 
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and 



Jt = 2.4 + O.67I/2 ^ 0.5I7 + 1.257 



3/2 



fi 



8.67^ + 1.357^/2 
225 - 177 + 72 ' 



(2.4) 



(2.5) 
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Whereas the plasmon decay depends in a comphcated way on the photon dispersion 
relation in media, its main features can be understood in an approximation where the 
photons are treated as particles with an effective mass equal to the plasma frequency, 
Up, which in the zero-temperature limit is given by |1| = Airane/ Ep, with a the fine- 
structure constant, the electron density and Ep the Fermi energy of the electrons. Up 
is of the order of a few tens of keV in WDs, slightly higher than the typical WD internal 
temperature, which is a few keV [l]. For the plasmon decay to happen, the decay products 
have to be kinematically accessible. Thus, when we talk about new light particles in this 
context we mean particles lighter than a few tens of ke V. 

It is not immediately obvious how the energy loss of Eq. (2.3) compares to the previous 
ones because of its complicated form, but the differences can be easily visualized in the 
simple plot in Fig. |4| For WDs whose internal temperature is below 4— 5x 10^ K, the cooling 
is dominated by photons, and perhaps axions. Above that temperature, the plasmon decay 
into two light particles becomes the main source of energy loss. The contribution from the 
bremsstrahlung of a pair of fermions is always negligible on the plot. 

It is useful to translate from temperature to Mboi. From Eq. (2.1), multiplying by 
a typical WD mass, M-wd? that we take to be 0.6 solar masses, we obtain the photon 
luminosity, = My/De^. Plugging it into Eq. (1.2) we obtain an expression that relates 
the temperature Tj to Mboi. Thus, a temperature of 4 — 5 x 10^ K corresponds to values 
of Mboi between 6 and 7, which is indeed where we see the neutrino dip in Fig. [T| The 
plasmon decay into a pair of light particles constitutes the dominant cooling mechanism 
for Mboi < 6 — 7, the exact figures depending on the mass of the WD and the properties 
of the envelope. 

Note that in models where a light boson couples to the electrons through a Yukawa cou- 
pling, the important cooling mechanism is the bremsstrahlung where the boson is produced 
singly, as in Fig. [s] (a). Such is the case for the DFSZ axion. In other models, instead, 
the light bosons couple indirectly to the electrons through a mediator and can only be 
produced in pairs, as for example when the bosons are charged under a new symmetry. 
This is the case for models with a dark sector, for instance, which we study in section |4} 
The dominant production for these bosons is then no longer the bremsstrahlung, but the 
plasmon decay. 
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3. A generic constraint on models with new light particles 



In this section we discuss generic constraints from WD cooling due to plasmon decay into 
new light particles, that can be either fermions or bosons, the only requirement being that 
they be lighter than a few tens of keV, in order for the decay to be kinematically possible. 
As mentioned at the end of the previous section, such a process affects the LF for values 
of Mboi below 6 — 7. Particle physics models in which new plasmon decay channels are 
open will potentially be in tension with the data, given the remarkable agreement between 
standard cooling mechanisms, that include neutrino emission, and the observed LF [s]. We 
want to quantify how much the plasmon decay rate can deviate from the standard one, 
considering the neutrinos as the only decay products. 

To achieve this goal, it is useful to introduce a unified formalism reminiscent of the 
Fermi interactions for fermions. In order to compare with standard plasmon decay into 
neutrinos, it is necessary to describe the relevant interaction between neutrinos v and 
electrons e. The interaction is given by 

X = -^H^(l-75)^](e7.e), (3.1) 

where the contribution from the effective neutral-current axial coupling constant Ca is 
negligible for our purpose and can be ignored [16j. From this Lagrangian one can compute 
the plasmon decay rate into two neutrinos. The result is fTsllMlfTT] 



_ CyG\ ZsTTg . , 

where a is the fine-structure constant, Zs is the plasmon wavefunction renormalization, vr^ 
is the effective plasmon mass which enters in the dispersion relation u'^ — k"^ = 7rs{u,k) 
for a plasmon with frequency u and wave vector k, and the subscript s = {T,L} denotes 
the plasmon polarizations (transverse and longitudinal, respectively). The explicit forms 
for ttt and ttl are involved. They can be found, for example, in Ref. [l]. We just point 
out for this discussion that vr^ is proportional to a, so that T^^s goes to zero if we turn 
off the electromagnetic interaction, as expected. With the standard plasmon decay rate 
into neutrinos, Eq. (3.2), the energy-loss rate per unit mass is given by Eq. (2.3) with 
C^G^ = CyGp, i.e. 

(3.3) 



plasmon plasmon 



C^G^—CyG-p 



and the contribution to the luminosity that appears in Eq. (1.1) is simply Lj^ = Mwd t^^^"'^. 
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Let us now turn to new neutrino-like cooling mechanisms for WDs. For BSM models 
with new light fermions ip, the relevant interactions are given by 



= -C^G^(V^7^V^)(e7^e), or = -QG^(7/;7V^)(e7^e) (3.4) 

which are the appropriate analogs of the four-fermion interaction. The quantities and 
have been described above. For new light bosons which must be produced in pairs 
[see Fig. [s] (d)], the interaction is 

= -2QG<^(^0t^^0)(e7^e), (3.5) 

^ y 

with (j)^ d ^(j) = (j)^{d'^(f)) — (9^0"'' )0. The corresponding quantities are the effective coupling 
constant C<^ and the dimensionful parameter G^, which is the analog of the Fermi constant. 
For both interactions the plasmon decay into two new light particles is 

r^^=^p^^, (3.6) 

where {Cx,Gx} are given by {C^,G^} for new light fermions or {C^,Gtf,} for new light 
bosons. The extra plasmon decay channel will lead to an extra energy- loss rate per unit 
mass as in Eq. (2.3) with C^G^ = CxGx, i-e. 



plasmon plasmon 



(3.7) 

-CxGx 



and an extra contribution to the total luminosity given by = M-^d e^^^^™""^, as in 



Eq. (1.1). In the following the relevant constants will be denoted simply by {Cx,Gx} 
both for new light fermions and bosons. 

As already mentioned, in order not to upset the excellent agreement between standard 
WD cooling mechanisms js], i.e. from photon emission and neutrino emission (relevant 
only for hotter WDs), and observational data, we postulate that plasmon decay into new 
light particles must not account for more than the plasmon decay into neutrinos. In the 
massless limit, both for new particles as well as neutrinos, this constraint can be stated 



simply as [see Eqs. (3.2) and (3.6)] 

CxGx ^ CyGp. (3.8) 

In other words, any new sufficiently light particles {i.e. which are effectively massless in 
WDs), that can be produced through plasmon decay in WDs and can escape from WDs, 
generate extra cooling. This extra cooling must be subdominant compared to standard 
plasmon decay into neutrinos. To validate this order-one constraint, it is now necessary 
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to properly quantify the agreement between the standard coohng mechanisms and obser- 
vational data. 

The standard cooling mechanisms relevant for WDs are photon cooling and plasmon 
decay into neutrinos. Since we are interested in constraining models which lead to extra 
neutrino-like cooling, we focus here only on the dataset of DeGennaro et al. which covers 
bolometric magnitudes between 5.5 < Mboi < 12.5. This range is well understood and 
clearly exhibits the neutrino dip for Mboi around 6 — 7 (see Fig. [s]). Moreover, the dataset 
of DeGennaro et al. has the smallest error bars in this range and only contains DA WDs. 

We start by minimizing the for the LF assuming standard cooling mechanisms and 
Mestel's approximation. The free parameter is the WD birthrate. The best fit implies 
a birthrate ~ 1.6 x 10^^ pc^'^ Gyr^^, which is a reasonable local WD formation rate. 



for x^in ~ 24.9. Since there are A^exp = 18 data points and A'th = 1 free parameters, 
this provides a decent fit with a reduced chi square Xredmin — Xmin/^dof = 1-47, where 
Adof = A'cxp — iVth = 17 is the total number of degrees of freedom. 

Next we determine the 90% confidence level exclusion contours for extra cooling from 
plasmon decay into new light particles, assuming the latter are massless. Since the new 
plasmon decay channels are reminiscent of the standard plasmon decay into neutrinos, we 
take here = SxL,y, where Sx determines the ratio of the new extra luminosity to 
the neutrino luminosity L^. Now we take Sx as our only free parameter, leaving the WD 
birthrate fixed to the value determined above. Thus we still have N^^of = 17. We then 
compute the new chi square, x^i including the Lx contribution. From Fig. 36.1 in Ref. [Is] 
we find that must be such that 

Ax' = - xLn < 24.8 (3.9) 

otherwise the extra cooling is excluded at 90% confidence level. Imposing the condition 
Ax^ < 24.8 translates into the constraint 

5. = ^= (^)'<0.99, (3.10) 



which is equivalent to the one in Eq. (3.8), obtained from a simpler and more intuitive 
physical argument. 

We provide in Fig. [5] three curves for the LF obtained from realistic models. The top 
one includes only standard cooling (with neutrinos), while the two lower ones include extra 
neutrino-like contributions with Sx = 0.5 and Sx = I respectively. 
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Fig. 5: Theoretical luminosity function for WDs. The curves shown include the different 
contributions to Eq. (1.1) and correspond to values of Sx = 0,0.5,1 from top to bottom 
for the = SxL^ contribution. They are superimposed on the data points by DeGennaro 
et al. i6j. 



4. Three examples 

In this section we consider three examples of BSM scenarios. The first two are super- 
symmetric extensions of the Standard Model (SM): in the first, the light particle is the 
neutralino, while in the second, it is the axino. We show that WDs do not put competitive 
bounds on these models. The situation is different in the third example, where we consider 
models with a dark sector, in which case the WDs bounds are very competitive. 



4.1. A light neutralino 

The neutralino xo is often the lightest supersymmetric particle in the Minimal Supersym- 
metric Standard Model. It can be very light, even massless, and still evade all current 
experimental constraints (l9]. For the production of light neutralinos in WDs, that would 
predominantly occur via plasmon decay, we consider the four-fermion interaction obtained 
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e" > 1^ > Xo, 1pa 

1 e 

e+ ^ * ^ Xo, 

(a) 

Fig. 6: Processes for the production of light neutralinos or axinos in WDs. The first diagram 
represents the relevant production mechanism for plasmon decay into neutralinos or axinos 
through a selectron exchange. The last diagram corresponds to the first diagram where 
the selectron is integrated out. The red dot is the corresponding four-fermion effective 
interaction. 




Xo, Wa 




Xo, 



(b) 



from integrating out the selectron e (see Fig. |6j) 

-^xo = -C^,G,{xoYfxo){e^,e\ (4.1) 



where Gg = -. — i-?, 1 and Cy,, = f 
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with e the electric charge, 6w the weak mixing 



angle and mg the selectron mass. Since Gp = „ . /f^^ 2 where mw is the W gauge boson 

O Sin (7 TTl 



mass the constraint Eq. (3.8) can be translated into a lower bound on the selectron mass 
of 



me > I ^^-^ 1 tanOwmw = 45 GeV, (4.2) 



where mw = 80.4 GeV and sin^ 9iy = 0.23. Thus, in order to have a significant impact 



on the LF one needs a selectron lighter than the W gauge boson. The bound in Eq. (4.2) 
applies to the case of a massless neutralino. Turning on a small neutralino mass has the 
effect of pushing the WD bound down to even lower selectron masses. Such light selectrons 
are already excluded by LEP searches. Note that supernovae, contrary to WDs, provide 
a better arena to constrain the mass of a light neutralino [2l]. Nevertheless, WD cooling 
bounds do not seem competitive for this process. 



4-2. A light axino 

A light axino is in principle very interesting in this context. It has already been argued 
that the inclusion of an axion gives a better fit to the LF. If supersymmetry (SUSY) is 
realized in nature, the axion would be necessarily accompanied by its fermionic partner, the 



very light neutralino, to^„ ^ 1 GeV, is almost purely bino and does not couple to the Zq 
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axino, which could also be very light (see e.g. [22]). The axino could be pair-produced in 
the plasmon decay and contribute to the high luminosity part of the LF. When combined 
with the contribution of the axion one might hope to get an even better fit. Unfortunately, 
as we explain in the rest of this section, the axino interacts way too weakly so that its 
contribution to the LF turns out to be completely negligible. 

Recall that the coupling of axions to electrons is given by igae'-f^e, with g = rrie/fpQ. 
In SUSY there is a corresponding axino-electron-selectron interaction that can be written 
as igeeipa, where ipa denotes the axino. If we integrate out the selectron (see Fig. |6]), 
the resulting four-fermion interaction between two electrons and two axinos is scalar-like 
[e.g. (V'aV'a)(ee)] instead of vector-like [e.g. (V'aT^^a) (67/^6)] and thus does not even allow 
plasmons to decay to pairs of axinos. Being more precise and starting from the deriva- 
tive interaction between the axion and electrons instead, one obtains higher-dimensional 
operators after supersymmetrizing and integrating out the selectron, i.e. four-fermion in- 
teractions between two electrons and two axinos with extra derivatives, which are thus 
temperature-suppressed compared to the usual plasmon decay. Most importantly how- 
ever, these interactions are always at least suppressed by g^, which is incredibly tiny for 



reasonable /pq ~ 10^ — 10^^ GeV. Therefore, although the constraint Eq. (3.8) cannot be 
applied directly here, the universal suppression just mentioned makes a possible production 
of axinos absolutely unobservable in WDs. 

4.3. A dark sector 

As seen in the two previous examples, WD cooling might not seem to lead to any strong 
bounds on new light fermions. The situation is however much more interesting when one 



considers models of BSM with massive dark photons 23 . In these models, which could be 
of relevance as models of dark matter, a dark sector communicates with the SM =^sm 
solely through kinetic mixing =SfsM(g)D l24j, i.e. 



— =S^SM + -^D + -^SM^D, where ^sm(^d — -^F'^F^ ■ (4-3) 

Above the electroweak scale the kinetic mixing occurs with strength sy between the hyper- 
charge gauge group f/(l)y, with the corresponding = d^B^ — duB^, and a new Abelian 
gauge group U{1)d, with F^^'^ = d'^A'^ — d'^A^^, where A'^ is the f/(l)D gauge boson, i.e. the 
dark photon. Below the electroweak scale the mixing involves instead the electromagnetic 
gauge group, and e = ey cos 9w- The dimensionless parameter e, which should be gener- 
ated by integrating out massive states charged under both SM and dark gauge groups, is 
naturally small, e ~ 10~^ — 10~^. Thus, after rotating the fields appropriately such that 
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(a) (b) 

Fig. 7: Process for the production of light dark sector particles in WDs. The first diagram 
represents the relevant production mechanism for plasmon decay into hght dark sector par- 
ticles through a dark photon exchange. The last diagram corresponds to the first diagram 
where the dark photon is integrated out. The red dot is the corresponding dimension 6 
operator. 



gauge bosons have canonically-normalized kinetic terms, the SM fields become millicharged 

i.e. 



under the dark gauge group 25 



^SM^D = -eeJ^'^A^, (4.4) 



where is the SM electromagnetic current. 



Thus, in models with massive dark photons, WD plasmons could decay, through off-shell 
massive dark photons, to light dark sector particles if they are kinematically available. Note 
that both dark photon decays to bosons and fermions result in two-particle final states. 
Thus such plasmon decay through massive dark photons into light dark sector particles is 
reminiscent of plasmon decay into fermions [e.g. neutrinos) irrespective of the spin of the 
light dark sector particles [see Fig. [s] (c,d)]. Therefore the relevant constraints for plasmon 
decay in models with massive dark photons are equivalent to the constraint discussed in 
section [3l 

In order to determine the resulting constrained parameter space it is necessary to inte- 
grate out the dark photon, as shown in Fig. [7) This leads to the interaction 

= -GcJ^r ^ -GBlC^ij^i^) + 2Q(^0t^A^0)](e^^e), (4.5) 

where the dark constant is Gd = and = Q^, = Here m^j, is the dark 

photon mass, is the dark fine-structure constant and Q^,^ are the dark particle charges 
under the dark gauge group. Note that dark photon decay into a pair of dark gauge bosons 
is generically not kinematically accessible because the masses of the dark photon and of 
the other dark gauge bosons are usually of the same order, as, for example, the Z and W 
gauge bosons in the SM. 
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10-1 



Fig. 8: Parameter space exclusion of dark forces with light hidden sector particles from energy 
losses in WDs. The blue shaded region is excluded by WD coohng for dark sectors with 
hght (< few keV) dark particles. For experiments, which usually assume the dark photon 
decay is predominantly into the SM, shaded regions correspond to completed direct searches 
while curves show future reach. For the electron and muon anomalous magnetic moments, 
shaded regions are excluded by measurements. The reader is referred to the text for more 
details. 



Comparing with plasmon decay to neutrinos as discussed in section |3j the constraint 



Eq. leads to 



2 < 



1.09 X 10 



-10 



rriAr 



V 



(4.6) 



IGn^a VGeV/ ' 

C^^tf). In Fig. [s] we show the constraint Eq. (4.6) and regions in parameter 

i.e. 

colliding experiments: 



where Cd 

space which have already been explored or will be explored by future experiments l26 
beam dump experiments at SLAC: E137, E141 and E774 27 ; e~^e 



BaBar 26 28 and KLOE |29); and fixed-target experiments: APEX 30, DarkLight 31 



HPS [32], MAMI f33| and VEPP-3 34 . Fig. |8| also shows excluded regions from electron 

(tte) and muon (a^ 

" 1111. . 1 ~ a, 



anomalous magnetic moment measurements 35 
For reasonable dark sector parameters where aj) 



one has Cd ~ 1 and 
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10~^ — 10~^, and thus all but a small fraction of the relevant dark sector parameter space 
is excluded by WD cooling and the above-mentioned experiments become obsolete if dark 
photons couple to new dark sector fermions and/or bosons which are effectively massless 
in WDs, i.e. lighter than a few keV. In other words, to be viable models of dark photons, 
any model probed by the above-mentioned experiments cannot have dark sector fermions 
and/or bosons lighter than a few keV due to WD cooling. 

Note however that all of the experiments shown in Fig. |8] — apart from DarkLight, 
VEPP-3 and the anomalous magnetic moment measurements — assume that dark photons 
decay predominantly back into the SM. Although this is not possible in WDs (dark photons 
could only decay back into electron-positron pairs which are not kinematically accessible), 
this assumption forbids either light dark sector particles, in which case the WD constraint 
presented here is irrelevant; or large dark fine-structure constant (relative to ae"^), for which 
dark photon decay rate into invisible channel dominates. 

To investigate this last possibility, we include in Fig. |8] (see dashed blue line) the WD 
cooling constraint for which the dark photon decay rate into visible channels dominates 
over the decay rate into invisible channels, i.e. Finvisibic ^ Tvisibie or C^od ^ It is 

interesting to see that the experiments which are sensitive to invisible dark photon decays, 
i.e. DarkLight and VEPP-3, are still tightly constrained by the WD cooling even when 
dark photons decay predominantly back into the SM. 

Note that the constraint must be modified for a very light dark photon (again lighter 
than a few keV), since it could be produced on-shell, which would result in an enhancement 
of the cooling rate. The resulting constraint would then be even tighter. For such a light 
dark photon bremsstrahlung might also become important. Finally it would be of interest 
to study astrophysical cooling constraints from more energetic objects, like supernovae, to 
relax the restriction on the masses of the dark particles produced. 

5. Discussion and conclusion 

We studied constraints from the WD LF on BSM models with new light particles. When- 
ever these light particles are produced in pairs, whether they are fermions or bosons, the 
dominant production mechanism in WDs is (usually) given by the plasmon decay. Such a 
decay is responsible also for the production of neutrino pairs, whose effect is well understood 
and clearly visible through the dip at Mboi ~ 6 — 7 in the LF curve. Adding a significant 
decay into new light particles would deepen the dip, which would then be in disagreement 
with the data. This constrains part of the parameter space of these BSM models. More 
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quantitatively, one needs to compare the strength of the interaction between the new hght 
particles and the electrons with the interaction between neutrinos and electrons, i.e. the 
Fermi constant Gp, and require that the former do not exceed the latter. 

We applied this constraint to three models. We first consider a supersymmetric model 
with a light neutralino and showed that the WD constraint is not competitive with existing 
collider bounds. The situation is analogous with an axino, whose interaction is even further 
suppressed with respect to the neutralino, and does not lead to any interesting constraint. 
We then explored models with a dark sector, for which the bounds are much more inter- 
esting. That is due mainly to the fact that the dark photon, that mediates the interaction 
between the electrons and the light dark sector particles, can be light (~ ^(GeV)), which 
enhances the plasmon decay rate. It turns out that the limits on the dark sector parameter 
space from energy losses in WDs, as shown in Fig. |8j are extremely competitive and render 
other experiments obsolete if the dark photon couples to light (~ ^(keV)) dark sector 
particles. Said differently, the dark photon models which are probed by these experiments 
cannot have light dark sector fermions and/or bosons, due to WD cooling. 

Note: During completion of this work An et al. [36^ posted a paper on stellar constraints 
for dark photons. There is no overlap between our work and theirs since they consider 
dark photons with hard Stiickelberg masses. 
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